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ABSTRACT 


When liquid is sealed at a high temperature, it flashes inside the seal due to the 
pressure drop and/or the viscous heat dissipation. Typical examples of appli- 
cations are light hydrocarbons in petroleum refineries, hot water in boiler feed 
pumps and reactor coolant pumps, and cryogenic fluids like liquid oxygen or hy- 
drogen in rocket turbopumps. The two-phase seals generally exhibit more erratic 
behavior than their single phase counterparts. Thermal effects, which are often 
neglected in single phase seal analyses, play an important role in determining 
a seal behavior under two-phase operation. It is necessary to consider the heat 
generation due to viscous shear, conduction into the seal rings and convection 
with the leakage flow. Analytical models developed by the past researchers work 
reasonably well at the two extremes - for low leakage rates when convection is 
neglected and for high leakage rates when conduction is neglected. None of the 
two models addresses the intermediate leakage rate problems when conduction 
and convection are both important. Most 6eals, however, operate in this range. 

A preliminary model, known as ‘Film Coefficient Model’, has been presented 
here which considers conduction and convection both, and allows continuous 
boiling over an extended region unlike the previous low-leakage rate model which 
neglects convection and always forces a discrete boiling interface. 

Another simplified, semi-analytical model, based on the assumption of isother- 
mal condition along the seal interface, has been developed for low leakage rates. 
This is basically a modification of the previous low leakage rate model and the 
range of applicablity of both is somewhat limited. The seal performance curves 
can be obtained quite easily. However, the ‘Film Coefficient Model’ may be used 
for more accurate and realistic description. 

Under certain circumstances, seals seem to exhibit self-sustained oscillations 
even when all the applied conditions remain steady. This was observed as ax- 
isymmetric fluctuations in the film thickness accompanied with periodic interface 
temperature variations. The two models, described earlier, will not be able to 
predict this anomalous seal behavior. A preliminary transient heat conduction 
model has been presented here which might help building a comprehensive model 
to describe the oscillatory seal behavior. 

The advantages and limitations of these models are discussed and some sug- 
gestions for future works are given. 
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Chapter 1 

INTRODUCTION 


Seals are mechanical devices used to restrict leakage of fluids, for example, when 
a rotating shaft penetrates a stationary housing which encloses the pressurized 
fluid. The tolerable leakage rate depends on the nature of the sealed fluid; leak- 
age of expensive, toxic, corrosive, explosive or flammable fluids must be reduced 
to a minimum. The life and reliabilty of seals are also of major concerns among 
the users to reduce equipment and process downtime. Somtimes when system re- 
dunduncy is kept at a bare minimum, for example in airborne and space vehicles, 
seal failure could cause serious system malfunctioning. 

In January 1986, the whole world suddenly became aware of the crucial im- 
portance of fluid sealing technology when the US shuttle “Challenger” tragically 
exploded shortly after leaving the launch pad. A joint sealed by rubber O-rings 
had failed. This episode had the characteristics of many a sealing problem. The 
component involved was of relatively low-value in its own merit, but the con- 
sequential cost of failure was totally disproportionate. The failed O-ring was 
a static seal and much less complicated in operation than the dynamic seals 
discussed here. 

The Figure 1.1 [1] gives an overview of different types of industrial sealing 
devices available. Among all kind of seals, the ‘Mechanical End Face Seals’ (also 
simply called ‘Face Seal’) are the dominating category of major industrial seals 
and have been given the special and extensive considerations. These are axial 
sealing devices with fluid leaking in the radial direction. The fluids need to be 
sealed, range from water, petrolium products, oil, natural gas, air and toxic chem- 
icals to cryogenic fluids like liquid oxygen and hydrogen (Space Shuttle Turbo 
Pumps). These seals handle pressure up to 5000 psi (~350 atm.), temperatures 
from up to 1000 °C down to -250 °C and a rotational speed of up to 60,000 RPM. 

A very good treatise on mechanical face seal designs, basic configurations, op- 
eration and lubrication mechanisms has been given by Ludwig & Griener [2,3]. 



Figure 1.2 shows the schematic diagram of a face seal. The primary sealing is 
accomplished by a nonrotating ring ( called primary seal ring or stator) that 
bears against the face of a rotating ring (called seal seat or rotor) mounted on 
the shaft. Occasionally co-rotating and counter-rotating seals (advanced aircraft 
engines) are encountered where both the rings are rotating. Between the stator 
and the housing, there are multiple springs which give it the flexibility in the ax- 
ial and two angular modes about orthogonal diametrical directions. Secondary 
seals are provided between the stator ring and the housing. Typically these seals 
are elastomeric O-rings. They self-energize under pressure and tend to fill in 
tlie asperitites and voids on the surfaces in contact and hence minimize leakage 
through secondary sealing surfaces. Sucessful operation of seals requires satis- 
faction of seemingly competing demands. In order to reduce wear and maintain 
integrity of the sealing surfaces, it is desirable, if not essential, to achieve and 
maintain seperation of faces by a lubricating film. At the same time face seper- 
ation must be kept extremely small (~ 2-3 /im) in order to minimize leakage. 
These requirements must be dynamically met in changing operating conditions 
and in the presence of machinery vibrations. 

Figure 1.3 shows the typical forces a stator experiences. This is an ‘outside- 
pressurized’ arrangement with the high pressure fluid at the Seal OD. This con- 
figuration offers a few advantages as the centrifugal forces tend to retard leakages 
and centrifuge solid particles away from the sealing surfaces giving a self-cleaning 
feature. The sealed fluid leaks through the gap between the seal rings and pres- 
sure drops due to friction and inertia. A typical pressure profile, P(r), for an 
axisymmetric gap is shown in the figure. If the gap is nonaxisymmetric, hydro- 
dynamic pressure will build up inside the seal due to relative tangential motion 
of the seal faces and the pressure profile becomes also a function of the circum- 
ferential location i.e P(r,$). For an axisymmetric gap, there is no hydrodynamic 
pressure generation and the pressure distribution inside the seal is same as if 
both the seal faces are stationary (if centrifugal inertia effects are neglected). 
The pressure distribution, so obtained, often referred to as hydrostatic compo- 
nent.. This fluid pressure P , tends to open the seal gap. On the other hand, 
the axial loading from the sealed fluid pressure and the spring force ,F, , acts 
behind the stator and tends to close the seal. The expressions for the opening 
and closing forces are given below. 

F 0 — J J P(r,6)rdrd9 

Pc = *(r 2 0 - r£,,)P 0 + F, + 7r(r^, - rj )P { 

‘ r boi' is called the ‘balance radius’ by changing which closing force ,F C , can be 
controlled. If the closing force, F c , is equal to the opening force, F 0 , at the 


operating point, then a lubricating fluid film is maintained at the seal interface. 
Under this situation, the seal operates in a ‘non-contacting mode’ and is called a 
‘balanced seal’. The corresponding clearance is called the ‘operating clearance or 
film thickness’. On the other hand, if the closing force is greater than the available 
opening force, the asperity contact takes place at the seal interface and the force 
balance is achieved with the help of the mechanical contact pressure. In this 
case, the seal operates in a ‘contacting mode’ and is called an ‘unbalanced seal’. 
The contacting seals are supposed to operate with a minimal contact pressure, 
otherwise heavy wear at the surfaces would cause prematured seal failure. These 
seals are generally used for low to moderate pressure services and noncontacting 
seals for high pressure applications. In chemical and petrochemical industries, 
contacting mode is primarily chosen for almost all sealing applications to reduce 
leakage of hazardous fluid as much as possible even at the expense of seal life. 
The third situtation arises when the opening force exceeds the applied closing 
force. In that case, the seal pops open causing high leakage and seal failure 
is said to have taken place. For a given design with a certain balance radius, 
rtaf, the closing force is constant for a given operating pressure whereas the 
opening force is dependent on the gap geometry and speed. The information 
most useful to the seal designers is the ‘Opening Force vs. Nominal Clearance’ 
curves, typically known as *F-h’ curves, for different speeds and system pressures. 
One necessary requirement for a stable and sucessful seal operation is to have a 
negatively sloped ‘F-h’ curve around the operating point (which means positive 
film stiffness), otherwise seal faces will collapse and give unacceptable contact 
load and rapid wear. Examples of typical ‘F-h’ curves are shown in Figures 1.9 
and 1.10 and they are expalined later. 

The seal lubricating film is usually very thin (in the range of few microns) 
and therefore, very small irregularities, thermal and pressure distortions, and 
face runout motions can have a dramatic effect on seal performance. Thus the 
primary seal cannot, in general, be visualized as two perfectly fiat and parallel 
surfaces. Some possible geometries are illustrated in Figure 1.4 [2]. The waviness 
(geometry a) and angular misalignment (geometry b) are most likely sources of 
hydrodynamic pressure bulid-up. Coning (geometry c) affects the hydrostatic 
pressure distribution and film stiffness. Externally imposed axial vibration (ge- 
ometry d) can produce squeeze film damping. Parallel misalignment (also called 
radial eccentricity) and shaft whirl (geometry e & f) impart a radial velocity 
component to the fluid particles which can affect the leakage. Some of these 
seal geometries, particularly the angular misalignment (geometry b) introduces 
dynamic forcing function with frequency same as the shaft rotation rate on the 
flexible ring which would consequently exhibit oscillations in axial and angular 
modes. One particular interest to the seal designers is to whether the flexible ring 



would be able to dynamically track the rotor without metal to metal contact, 
for a given amount of rotor misalignment (commonly called ‘runout’). In the 
‘dynamic tracking’ analysis, the fluid film is modeled as nonlinear springs and 
dampers. The damping comes from the ‘sqeeze film’ effects. The elastomeric 
O-ring also provides considerable amount of stiffness and damping both of which 
are frequency dependent. If O-ring starts slipping, Coulomb damping also comes 
into play. It has been found that higher the film stiffness is, better the dynamic 
tracking ability of the seal ring becomes. 

The engineering of seal can involve fluid mechanics, heat transfer, elasticity, 
thermodynamics - equilibrium and nonequilibrium, statistical mechanics, dy- 
namics, chemistry and metallurgy, to name a few of the most frequent areas of 
concern. Usually each effect can be analyzed by itself, but then the integrated 
effects must be evaluated for a complete analysis of a sealing system. As indicated 
previously, seals are characterized by surfaces in relative motion seperated by a 
very narrow gap. In order to ensure proper operation, very small differences in 
the dimensions of the seal part must be maintained while in operation. Deforma- 
tions in geometry due to imposed thermal gradients, frictional heating, pressure 
and mechanical contact forces must be held to a minimum. In any cases, net 
deformation must be no more than microvalues. Depending on the imposed con- 
ditions, seals operate basically in three different lubrication flow regimes shown 
in Figure 1.5 [2]. The ‘non-contacting’ seals usually operate with ‘full film lubri- 
cation’ whereas the ‘contacting’ seals can operate either in the ‘boundary (also 
called ‘mixed friction’) lubrication’ or ‘dry sliding lubrication’ regime depending 
on the excess maginitude of the closing force over the available opening force. 

If the sealed fluid is a gas (usually operating in a noncontacting mode), some- 
times the ‘mean free path’ of the molecules may be of the same order or more 
than the nominal operating seal clearance (Knudsen No. > 1) in which case 
continuum fluid mechanics with no-slip boundary conditions is no longer valid 
and ‘slip flow theory’ and ‘statistical mechanics’ are needed to describe the fluid 
flow. 

The selection of materials for seal rings is also a very important aspect of 
seal design. It requires extensive tribological testing to come up with a suit- 
able material combinations for a certain specific application. In general the seal 
ring materials should have good mechanical and thermal shock resistance, wear 
characteristics, corrosion resistance, self lubrication property, high modulus of 
elasticity, to name just the important ones. Carbon-graphite usually meets most 
of the requriements. It is quite frequently chosen in combination with some other 
compatible hard material like tungsten or silicon carbide. There are other combi- 
nations of seal materials used also e.g carbon-graphite vs. stainless steel, tungsten 
carbide vs. tungsten carbide etc., depending on the operating conditions. 



With all the complexities and highly coupled effects that govern a seal be- 
havior, it is not a secret that reliable and accurate design analysis for face seals 
does not exist. Any new seal design must be tested in the laboratory because 
a prediction of eventual performance is not possible on purely theoritical basis. 
Again there is a wide variation in performance of seals of “identical design”; a 
particular seal may fail after few hours whereas another seal belonging to the 
same class can last for several years or so. It is because of this reason seals are 
often termed as the ‘most unpredictable machine element’ used in the industry. 
However a good design and analysis tool is quite useful in evaluating one design 
against others. This procedure eliminates the need for building expensive proto- 
types and running time intensive laboratory tests for those designs which seem 
not so viable at the analysis phase [4]. Also modeling and analysis give more 
insight into the complex mechanism of seal behavior. Hence there have been 
quite a good amount of analytical and experimental work done in the face seal 
area for the last 25 years and efforts are constantly being made by the engineers 
and scientists to come up with better theoritical models for seal operation. Since 
it is extremely difficult to perform a comprehensive seal analysis including all 
different coupled effects, most of the analytical works are focussed on one or two 
aspects. As of yet some of the individual effects (e.g two phase flow modeling, 
high Knudsen number flow, mixed friction regime, wear model, nonequilibrium 
effects, to name a few) are not fully understood. 


1.1 Liquid Seals 

A number of investigators have analyzed face seal operating with incompressible 
liquid for different flow geometries. Etsion have extensively studied the angular 
misalignment effects on seal performance and stability. A misaligned face seal 
i6 shown schematically in Figure 1.6 [11]. He obtained a complete system of 
forces and moments acting on the flexible ring for different values of angular 
misalignment. These can then be used in a seal dynamic tracking analysis. 
Etsion [5] observed that any angular misalignment produces a radial force on the 
flexible ring which in turn causes a radial eccentricity. When this eccentricity is 
large enough, the pumping of fluid may take place which will affect the leakage. 
The seal coning however tends to reduce the magnitude of the radial force [6]. 
When the pumping takes place in a direction opposite to the hydrostatic pressure 
drop, it is known as ‘inward pumping’. This phenomenon was studied both 
analytically and experimentally by Findlay [7,8]. Analysis also showed that a 
flat outside pressurized seal with angular misalignment has negative axial and 
angular stiffness [9]. Also the hydrodynamic forces creates a transverse moment 



which leads the tilting moment by 90 degrees [10] that can cause seal wobble. 
However with coning, the stiffnesses might change sign depending on the relation 
between the angle of tilt and angle of coning [11,12]. For noncavitating flow, 
the effect of coning reduces the hydrodynamic transverse moment which would 
improve seal stability. The ‘narrow seal approximation’ is usually made in seal 
analysis for simplification. With this approximation, circumferential pressure 
gradient and seal curvature can be neglected. Etsion [13] compared the accurate 
results from numerical solutions with the approximate results and found that 
over a radius ratio, r,/r 0 , greater than 0.8, an accuaracy of less than 1% can be 
obtained and hence most cases this approximation is justified. 

For low pressure and/or high speed seals, lubricant cavitation is possible due 
to hydrodynamic effects. This has been experimentally observed. An interesting 
work on this subject has been published by Findlay [14]. The lubricant cavi- 
tation helps generating extra opening force because it prevents the generation 
of hydrodynamic pressure below the local vapor pressure of the liquid while not 
restricting the upper bound of the pressure. If cavitation did not occur, the com- 
ponents of hydrodynamic force would usually balance out and no net increase 
over the hydrostatic force would exist, which is not the case for low pressure 
seals. 

Sneck was one of the early investigators who made a very important contri- 
bution in the face seal analysis under incompressible flow. He published a series 
of papers [15] through [20] in ‘68 - ‘69 in which he addressed different aspects af- 
fecting seal performances e.g angular misalignment, radial eccentricity, tangential 
waviness, flow turbulence, centrifugal inertia and thermal effects. The centrifu- 
gal intertia term is included in the misalignment analysis in [15]. The centrifugal 
effects are shown to play a significant role in seal performance at higher speeds. 
For an outside pressurized seal, the regions of flow field may exist with a radially 
inward flow along the stationary surface and outward along the rotating one and 
under certain circumstances, there can be net zero leakage. The existence of such 
a region is a direct consequence of centrifugal inertial effects. This reverse flow 
phenomenon has been studied in detail in another paper [16]. The combined 
effects of misalignment and radial eccentricity is presented in [19]. The resulting 
leakage component can be inward (opposite to the direction of hydrostatic pres- 
sure drop) or outward (in the same direction as the hydrostatic pressure drop) 
depending on the phase angle between the misalignment and radial eccentricity. 
Sneck also studied eccentricity combined with surface waviness [20]. Again the 
direction of leakage component is shown to be dependent on the phase angle. The 
once per revolution waviness is found to be the main contributor in the pump- 
ing effect. Turbulent flow is analyzed in [17]. The turbulent nature of the flow 
is described by an isotropic apparent viscosity model and a power law velocity 




profile. The misalignment and surface waviness are found to be somewhat less 
influential with turbulent flow than with laminar flow. 

The analysis of face seal is often based on the isothermal flow assumption 
within the seal clearance. The validity of this assumption is usually argued on 
the basis that seal faces are often good conductors and hence will not permit large 
radial temperature variations. But even when the seal operates approximately 
isothermally, the temperature within the seal clearance need not necessarily be 
same as the cavity fluid temperature. An accurate prediction of seal performance 
requires an accurate evaluation of the fluid viscosity within the clearance space. 
A general thermal analysis procedure is presented in [18] to estimate the fluid 
operating temperature level inside the seal. No attempt has been made here to 
model the heat conduction through the seal rings. The upper and lower bound 
on the operating temperature are obtained by assuming adiabatic wall condition 
and zero thermal convection by the fluid respectively. In a recent review paper by 
Khonsari [21], an extensive survey of literatures pertaining to thermal effects in 
slider and thrust bearings is presented with summary of important contributors 
of leading researchers and designers. Since thrust bearings and seals have some 
similarity, this paper is referred here. One very common assumptions made by 
the seal analysts is to neglect the fluid temperature variations across the film. 
But viscosity variation across the lubricant film has been somtimes found to be 
responsible for generation of an appreciable load. In [21] many papers are cited 
which indicate the importance of transverse viscosity variation. King and Lauer 
[22] presented an experimental method by infrared spectroscopy to verify the 
existence of the temperature gradients through the film. 

Pinkus and Lund [23] also considered the effects of centrifugal forces in high 
speed seals. They mentioned that at the upper limits of laminar conditions, 
centrifugal forces reduce the load capacity considerably and alter the pattern of 
the lubricant flow. Koga & Fujita [24] included both the radial and centrifugal 
inertia terms in their analysis of high pressure water pump seals. They obtained 
better correlation of the analytical predictions and experimental results when 
inertia effects are considered than their previous analysis neglecting these effects. 

As mentioned before, the total closing force is supported by hydrostatic and 
hydrodynamic fluid pressures and often by partial contact of fluid faces (in con- 
tacting mode of operation only). For moderate to high pressure applications, the 
hydrostatic force component is predominant over the hydrodynamic component 
[25]. Since the film thickness is usually very small (of the order of few microns), 
any local surface deformations due to the interfacial pressure and the angular 
twist of the seal rings under pressure strongly influence the hydrostatic load sup- 
port and hence the seal performance [26]. For carbon rings with a relatively low 
modulus of elasticity, the distortions can easily be of the same order of magni- 
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tude as the nominal clearance of the seal. Thermal distorotions can also occur 
due to both axial and radial temperature gradients in the seal rings caused by 
the frictional heat generated at the interface. Any radial taper in the direction 
of the flow changes the hydrostatic pressure distribution and the film stiffness. 
A diverging seals (in the leakage direction) exhibits a negative axial stiffness 
which may lead to seal collapse and high wear. With the availability of fast 
computers and necessary softwares, the finite element (FE) analysis is commonly 
used for accurate predictions of pressure deformation and thermal distortion. A 
very important series of papers [27,28,29, to mention a few of them] have been 
published over the years by Metcalfe and his research staff at the Atomic En- 
ergy of Canada Ltd. (AECL) describing these analytical techniques. Analysis of 
seal-ring deflections due to applied pressure loadings, to thermal effects and to 
coloumb friction between components is described in [28]. The deflection sen- 
sitivity of seal components are expressed as ‘inflence coefficients’, evaluted with 
finite element analysis. He noted that coloumb friction gives rise to undesirable 
performance hysterisis when operating conditions are changed. The correlation 
of experimental results and theoritical predictions is presented in [29]. Salant 
[30] presented an analytical model of a generalized mechanical seal incorporat- 
ing the fluid dynamics of the film and the mechanical and thermal distortions 
of various seal components. He utilized the concept of ‘influence coefficients’ 
used in [28]. He found that the hydrodynamic forces due to waviness, roughness, 
misalignment and eccentricity produce insignificant opening force effects in com- 
parison with the available closing force for high pressure seals. The hydrostatic 
pressure is responsipble in carrying the most of the applied load. Hence the face 
deformation, particularly ‘coning’ is the most likely controlling mechanism for 
load support for these kind of seeds. Based on this idea, Salant presented a novel 
design of an electronically controlled seal in [31]. A microcomputer based real 
time control system and electro-mechanical actuator dynamically adjust the seal 
coning and hence the film thickness, based on information received from the sta- 
tor which monitors conditions of the film. This arrangement can greatly reduce 
face contact while limiting leakage by continuously optimizing film thickness. 
This would lead to a reduction in seal damage and wear and increase in seal life. 
Li [32] presented a finite difference heat conduction model for calculating the 
temperature distribution in the seal rings and resulting deformations. He con- 
sidered one dimensional toroidal deformation model in which the seal rings will 
preserve the geometry of the radial crosssection after deformations. Doust and 
Parmar [33,34] numerically analyzed axisymmetric distortions due to pressure 
and thermal effects using ‘boundary integral element’ (BIE) method and corre- 
lated the results with the experimental measurements. They remarked that the 
BIE method is substantially more economical in terms of both computing time 
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and storage than FE for the same lavel of accuracy. The main object of their test 
was to measure the fluid film geometry using capacitance type proximity probe, 
as a function of pressure. The sealant pressure and thermal effects essentially 
caused toroidal rotation of the faces for those seals used by them. This rotation 
rate was found to be fairly insensitive to the interface pressure profile. They 
also observed hysterisis effect due to secondary 6eal friction. In a recent paper 
[35] by them, the effects of thermoelastic transients have been presented. The 
transient thermal distortion can be an order of magnitude greater than that at 
steady state. Transient response is worse for a shorter section than a longer one, 
although the time to reach steady state can be more than an hour for a long 
seal component. This is an interesting work since field surveys do suggest that 
transient operations can be more detrimental to the seal life than steady state 
running. 

There are not as much work done in the ‘mixed-friction lubrication’ area for 
the contacting seals as in the ‘full-film lubrication’ regime for the noncontact- 
ing mode of operation, mentioned so far. The obstacle to further advancement 
in contacting seal technology is that relationship between controllable design 
parameters and performance parameters are not well understood. Lebeck has 
published a number of papers [36,37,38,91] on ‘mixed-friction’ flow modeling and 
contacting seal analysis. He developed a model [36] which takes into account load 
sharing between mechanical and fluid hydrostatic pressure. The effect of wear 
is also modeled in order to predict how the radial profile alters and influences 
the hydrostatic pressure distribution with time. The experimental evaluation of 
the model i6 reported in [37]. In a contacting seal operation, an unstable phe- 
nomenon is observed by a number of researchers, including Kenndy and Grim 
[39], in which case a very slight amount of initial waviness on seal faces grows 
during seal operation. When this unstable condition, called ‘thermoelastic in- 
stability’ occurs in an operating face seal, the consequences • nonuniform wear, 
accentuated waviness, and high localized stresses and temperatures - can be very 
detrimental to seal performances. Kiryu et al. [40,41] reported the generation of 
a “ringing” sound in a contacting water pump seal. They attributed this phe- 
nomenon to self-excited vibration due to ‘stick-slip’ action, caused by transferring 
from fluid lubrication to dry sliding condition. Vibration mode in ringing sound 
generation is found to be mainly caused by the torsional and axial vibrations of 
the rotating shaft system. 

The preivious investigations, mentioned so far, are mainly based on steady 
state analyses. However, the angular misalignment is inevitably present on the 
rotating ring which introduces a dynamic forcing function on the flexibly mounted 
stator. Hence the ability of the stator ring to track the rotor in a controlled man- 
ner is of great importance for safe seal operation and as the demand for higher 
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operating speeds in rotating machinary increases, the importance of seal dynam- 
ics becomes more and more evident. A several researchers namely Etsion, Green, 
Metcalfe and others have addressed this issue analytically and experimentally in 
[42] through [55] A review of face seal dynamics covering the literature until 1981 
is presented in [50]. 

The flexibly mounted stator has basically three major degrees of freedom - 
axial and angular about any two orthogonal diameters. The twisting motion 
about axial direction is prevented by antirotation locks. If the radial stiffness of 
the O-ring secondary seal is low, which is usually not the case, then the stator 
can also move in the two perpendicular radial directions. The rotor transmits 
its angular motion to the stator via the thin fluid film seperating the two seal 
rings. For a given forcing function, the response of the stator depends on its own 
inertia, the stiffness and damping of the fluid film and the elastomeric O-ring. 
The fluid film damping comes from the squeeze effects. In some cases the squeeze 
effects are an order of magnitude higher than the combined hydrodynamic and 
hydrostatic effects and hence play an important role in dynamic behavior of 
face seals. The stiffness and damping coefficients of the fluid film ,both direct 
and cros6-coupled, in the three major d.o.f are calculated in [44,51] based on 
small perturbation theory. It has been found that narrower the seal, the less 
is the damping coefficients and at very small tilt, translational and rotational 
direct damping coefficients are an order of magnitude higher than the cross- 
coupled ones. The damping and stiffness characteristics of elastomeric O-ring are 
dependent on the amplitude and freqency of excitation and amount of sqeeze. 
The experimental determination of the O-ring dynamic properties are presented 
in [54,56,57]. With a large rotor misalignment, the stator response is usually 
large and sometimes sliding takes place at the O-ring interfaces and then the 
Coloumb friction becomes important. 

Dynamic analysis in [45,46,48,52] based on linearized small perturbation the- 
ory revealed three modes of operation : a stable mode in which a misaligned 
rotor is synchronously tracked by the flexibly mounted stator; a transition mode 
in which half-frequency wobble of the stator is superimposed on the previous 
synchronous tracking mode; and an unstable mode characterized by uncontrolled 
vibration of the stator, eventually causing failure. In the unstable mode, a seal 
will fail even with zero rotor runout. For low and moderate speeds, the stable 
mode seems to predominate. The stator tilt, however, differs from that of the 
rotor both in magnitude and direction. The difference and phase shift between 
two tilts result in relative angular misalignment between the rotor and stator. If 
this relative misalignment becomes too large, seal failure due to excessive leakage 
or even rubbing contact can occur even though the seal is dynamically stable. In 
[55], the complete nonlinear equations of motion of the stator are solved numeri- 


cally. The assembly tolerances in the form of initial stator misalignment and the 
dynamic properties of the elastomeric O-ring are accounted for in the analysis. 
Both statbility threshold and steady-state response of the stator are investgated. 
In general, it was found that the critical shaft speed corresponding to stability 
threshold is quite high. Hence, the dynamic stability should not be a problem 
in the majority of noncontacting seals. A more practical problem is the steady- 
state dynamic response of the stator resulting from rotor runout and assembly 
tolerances. The results of the numerical analysis were compared with those of 
the previous small perturbation analysis that provides much simpler closed form 
analytical solution. Very good correlation was found between the two analyses 
for most cases of practical applications. 

Etsion and Burton [43] observed self-exicted oscillations of seal ring in the 
form of precession and nutation. The wobble frequency was measured to be 
about 43% of the rotational frequency. Metcalfe [49] analyzed and tested a well- 
aligned face seal. He found that if the balance ratio is below a certain critical 
value the seal becomes hydrostatically unstable. If the elastomer stiffness in the 
tilting mode is insufficient to overcome this hydrostatic instability, the stator will 
exhibit wobble motion. The precession rate is theoritically found to be half the 
shaft speed if elastomer damping is insignificant (pure “whirl”) and progressively 
slower as damping increases. Etsion presented an experimental observation of 
the dynamic behavior of face seals in [53]. The forced response of the stator 
due to the rotor runout was monitored by means of three proximity probes. It 
was found that both the stator tilt and its phase shift with respect to rotor tilt 
are time dependent and vary synchronously with the rotor rotation. The time 
variation is attributed to the presence of two components of stator tilt. One 
component is fixed in magnitude but tracks the rotor tilt. The other component 
is fixed both in magnitude and direction and is due to nonaxisymmetric effects in 
the flexible support of the stator. As a result, the relative misalignment between 
the stator and rotor was found to be time dependent. The dynamically unstable 
seal behavior was also observed. At low supply pressure which means low film 
stiffness, results showed a sinusoidal perturbation at double the shaft frequency 
superimposed on the initial wobble due to angular misalignment. In the previous 
experiments, perturbation was always observed at half rather than double the 
shaft rotation rate. As the supply pressure was increased, this double frequency 
stable to unstable transition became a half-frequency transition instead. This 
higher frequency instability was not fully understood. 



1.2 Gas Seals 


The efforts on gas seal development started little later than the liquid seals and 
the works in this area are not equally numerous. In earlier time, the machineries 
having gases as working fluids, like gas compressors, used and some of them still 
use liquid seals with an oil-buffered arrangements. The reason for this that proper 
technology was not available to insure non-contacting mode of operation, which 
is absolutely essential for gas seals because of the poor lubrication properties 
and high speed of operation. The oil is kept at a pressure little higher than the 
sealed gas to ensure that only oil leakage could take place into the gas and seal 
would never run dry. In addition oil also leaks to atmosphere through another 
seal. Apart from the cost factor (about two orders of magnitude higher than 
the corresponding single phase seals), this design has some major disadvantages 
in terms of auxiliary equipments and space requirements. Also since product 
contamination with just a small amount of buffer fluid may create enough prob- 
lems, contacting liquid seals are typically used, which have inherently low and 
unpredictable life. Hence the need for a noncontacting seal development did 
arise for sealing gaseous fiiuds. Although the basic concepts are the same, the 
main difference between the liquid and gas seal analyses is that the governing 
equations describing the gas flow are nonlinear because of compressibility effects 
and inclusion of flow inertia, which are too important to ignore sometimes. The 
flow is often turbulent and choking may occur at the outlet. Also under the con- 
ditions of high velocity, the entrance loss effects cannot be neglected. The gas 
seals should also operate with a high film stiffness in order to have good dynamic 
tracking ability to prevent contact. As mentioned before, the pressure within 
the fluid film are generated hydrodynamically by the relative motion between 
uneven sealing surfaces and hydrostatically by frictional pressure drop through 
the seal. The hydrodynamic action ceases when the motion stops. Also there 
is no hydrodynamic pressure generation with parallel faces. To a limited extent 
all seal possess some hydrodynamic characteristics as a consequence of geometric 
imperfections and unplanned unevenness such as inherent or pressure induced 
circumferential waviness or micro-irregularities. These effects are usually quite 
small. The hydrostatic effects alone impart zero stiffness to a seal unless there 
is a radial coning in the flow direction. Because of these facts, some consious 
efforts have been made to enhance the hydrodynamic action rather than rely 
on chance variation, by having planned uneven hydrodynamic patterns on the 
seal surfaces. Some of the commonly used patterns are spiral groove, Rayleigh- 
step pads, radial grooves, as shown in Figure 1.7 [58]. These are called ‘hybrid’ 
seals. The hydronamic pattern is followed by a seal dam which offers restrictions 
to the fluid flow and most of the pressure drop takes place there. Because of 



hydrodynamic action, there are some areas of higher pressure and other areas 
with lower pressure. Figure 1.8 [58] shows the elevated pressure areas on the two 
seals. Figure 1 .9 shows the components of *F-h* curves for the hydrodynamic and 
the hydrostatic sections of a 'hybrid 1 seal with Rayleigh-step pads, analytically 
obtained by Shapiro [59]. The two curves must be combined to get the net film 
characteristics for the seal under consideration. No angular misalignment effects 
is considered in this analysis. It is evident from this figure that the hydrodynamic 
action indeed imparts a very high film stiffness, partcularly at small clearances 
and prevents the seal faces from touching each other. The closing force is usu- 
ally chosen so that the seal operates near the high stiffness region. Also it is 
seen that hydrostatic stiffness is almost zero and hence it does not contribute to 
the seal stability, although it may carry a major part of the closing force. Ex- 
periments performed by Ludwig [60] showed that seals with hydrodynamic pads 
outperformed the conventional seals used in small gas turbines. 

Some of the important researches on gas seals are documented in [59] through 
[72]. Cheng [61] analyzed few different designs and found that the spiral groove 
design gives higher stiffness than Rayleigh-step pad one. This same conclusion 
was also drawn by Sedy [58]. He also brought out an interesting point. As men- 
tioned before, most of the gas expansion takes place over the dam. The cooling 
effect associated with the expansion is sometimes several times more than the 
heating effect due to viscous dissipation. The net effect is the cooling of the gas 
near the dam and consequently a considerable amount of heat conduction takes 
place from the seal rings to the gas in the vicinity of the dam. The tempera- 
ture gradient, thus set up in seal rings, tends to distort the seal face in a way 
to produce a divergent flow passage which has an unstable effect and sometime 
causes seal contact at the outer diameter. Sedy suggested a wider dam design to 
overcome this problem because a wider dam would cause a higher heat genera- 
tion which in effect tends to neutralize the cooling effect due to gas expansion. 
Zuk [64] presented a quasi one-dimensional analysis for the flow of gas through 
seals. This model includes fluid inertia and entrance loses, in addition to viscous 
friction which is accounted for by a friction factor. Subsonic and choked flow con- 
ditions have been predicted and analyzed. This model is valid for both laminar 
and turbulent flows. Hsing and Carraro [74] used an efficient algorithim based 
on fourth order Runge-Kutta with adaptive step size to solve the same govern- 
ing differential equation. Shapiro [59] performed both steady state and dynamic 
analyses of a gas seal for jet engines. The seal dynamic response was found as 
a function of rotor misalignment and secondary seal friction. He theoritically 
obtained superharmonic response, about four times the shaft RPM, which has 
been confirmed experimentally. This phenomenon has been attributed to nonlin- 
ear characteristics of the O-ring. The flexible ring is found to loose its tracking 


ability if the rotor runout or the friction force is too large. 

If the gas inside a seal is at sufficiently low pressure, the molecular mean 
free path can become comparable to the film thickness. The fluid subjected to 
this condition does not behave entirely as a continuum fluid but rather exhibits 
some characteristics of molecular chaos. One may also expect to encounter these 
effects in regions having very sharp gradients of fluid properties such that these 
properties change appreciably in the space of a few mean free paths, regardless 
of whether or not the absolute density of the gas flow is especially low. The 
dimensionless ration, X/h (Knudsen number), is a measure of the degree of rar- 
efaction. When this ratio is large, the flow phenomena are mostly dictated by the 
molecular-surface interaction. This class of fluid flow is defined as “free- molecular 
flow”. For flows in which the value of Knudsen number is small, typically 0.01 
< X/h ~ 0.1, but not negligible as those in continuum mechanics, some depar- 
tures from the usual continuum flow phenomena may be expected to occur. The 
layer of gas immediately adjacent to the solid surface no longer assumes the same 
kinematic condition as the solid boundaries but has a finite relative “slip veloc- 
ity” and hence produces an apparent diminution in fluid viscosity. This is called 
the “slip-flow regime”. 

The few research works covering these non-continuum effects on lubrication 
are documented in [73] through [78]. Hsing and Malanoski [74] found that if 
lubricant is one of the gases having a large molecular mean free path, such as 
Helium, Neon or Hydrogen, the slip-flow phenomena could contribute substantial 
reduction in the performance of a thrust bearings, which is quite similar to face 
seals. Cans [75] derived a slip flow lubrication equation for an arbitrary Knudsen 
number from kinetic theory. Fukui and Kaneko [76] developed a more accurate 
generalized lubrication equation based on linearized Boltzmann’s equation. The 
experimental results obtained by them [77] agreed well with their numerical re- 
sults. Kubo et al. presented a finite element solution of the Boltzmann’s equation 
in [78]. 


1.3 Two-phase Seals 

When liquid is sealed at temperature higher than its saturation temperature at 
the outlet pressure, it flashes inside the seal due to the pressure drop and/or the 
viscous heat dissipation. Typical examples of applications where such two-phase 
flow may be encountered are light hydrocarbons in petroleum refineries, hot 
water in boiler feed pumps and reactor coolant pumps, and cryogenic fluids like 
liquid oxygen and hydrogen (LOX / LH2) in rocket turbopumps. The two-phase 
seals generally exhibit more erratic behavior than their single phase counterparts. 


The seals also have more stringent requirements in their performances because 
of severity in applications. As for example, light hydrcarbons are potentially 
flammable and explosive and hence certainly dangerous if allowed to leak. Since 
these hydrocarbons in gaseous phase are heavier than air, they usually form a 
thick dense cloud on the ground around the source. It constitutes a severe hazard 
[79]. In LOX / LH2 turbopumps, any seal failure due to excessive leakage can be, 
needless to say, extremely dangerous. Actually the face seals in the turbopumps 
failed repeatedly on the NASA test pads until they have been replaced with 
annular seals. Although annular seals are safe in operation, they allow very high 
leakage and did cost NASA millions of dollars worth of fuels during each testing 
and launching operation. At a later date, the face seals have been adopted 
sucessfully in the LOX / LH2 turbopumps for the Japanese H-l rocket [80]. 
Two-phase seal operation is also encountered in boiler feed pumps. It has been 
estimated that the boiler feed pump outages alone cost power companies several 
hundred million dollars each year in lost power revenues. It is believed that a 
high percentage of these problems is related to seal failures [97]. The reactor 
coolant pump (RCP) seals can also experience change of phase of the sealed fluid 
during station blackout conditions and exhibit excessive leakage. Failure of these 
precision components may result in a small loss of coolant accident (LOCA) in 
nuclear reactors [82]. 

Because of the severity of applications, the two- phase operation mechanisms 
are to be better understood in order to come up with a suitable design. The 
research works done in this area are reported in [79] through [102]. An inter- 
esting earlier papers on this subject was published by Orcutt [83]. He used a 
quartz runner to permit visual observation of the seal interface during operation. 
The experimental observations indicated the existence of a multiple phase film, 
charcterized by two large scale regions. The first region adjacent to the seal cav- 
ity was occupied almost entirely by water. The second annular region extends 
from the atmospheric edge of the interface to a semi-stable boundary with the 
liquid-filled region. This region was occupied by a mixture of liquid and vapor. 
The boundary moved towards the edge adjacent to the seal cavity with the rise 
in liquid and seal surface temperatures. Unstable operation was encountered 
with visible leakage as the cavity fluid temperature was increased. More than 
a decade later, Harrison and Watkins [84] and Wallace [79] reported a similar 
unstable two-phase operation with light petroleum products at elevated temper- 
ature. Under the unstable operation, the fluid Aim periodically broke down and 
reformed with violent fluctuations in torque. Seals showed both audible (while in 
operation) and visible (when taken apart) signs of distress. Seal operation was, 
however, stable at lower temperature. Barnard and Weir [85] reported seals oper- 
ating sucessfully with no visible leakage because of vaporization. The seal faces, 


they examined, all exhibited three concentric bands across their surfaces. Will 
[86] also observed the similar three banded appearance on 6ucessfully operating 
two-phase seals. No convincing causes are known. 

In theoritical modeling of two-phase flow through seals, most of the works 
so far have been done by Hughes and his coworkers and are reported in [87] 
through [93]. Basically two different models have been presented in [87,88,89] for 
low and high leakage flows respectively. The low leakage flow model is based on 
laminar flow and it considers heat conduction into the seal rings but neglects heat 
convection by the fluid and heat of vaporization. Boiling is assumed to be taking 
place at a discrete interface. The high leakage flow model is based on turbulent 
flow. This model disregards heat conduction through seal rings but takes into 
account convection, heat of vaporization, radial and centrifugal inertias. This 
model could predict continuous boiling over a finite region and also choking at the 
outlet under certain conditions. Beatty and Hughes [90] refined the turbulent flow 
model with better treatment of inlet losses. They obtained an anomalous ‘all- 
liquid choking’ situation in which the flow is choked but remains liquid all the way 
up to the seal exit. Beyond the exit, the liquid flashes immediately into vapor. 
Lebeck presented a mixed-friction model with phase change in [91]. He modified 
the flow equation for roughness effects and considered the load support due to 
asperity contacts. Hughes and his coworkers assumed an idealized semi-infinite 
heat conduction model whereas Lebeck used a more realistic seal geometry and 
boundary conditions and implemented a finite difference scheme to solve for 
seal ring face temperature distributions. Lebeck ’s model is valid for low leakage 
rates only. Beeler and Hughes [98] performed a dynamic analysis in the axial 
mode. They used the quasistatic ‘F-h’ curve obtained by using the adiabatic 
model to represent the fluid film. Squeeze film effects were ignored. With this 
limited dynamic model, they predicted self-sustained oscillations under certain 
conditions whereas failure due to metal-to-metal contact under other situations. 
Zuber and Dougherty [94] modeled the process of condensation and evaporation 
and derived a generalized lubrication equation. The two-phase region is treated 
as a dispersed homogeneous flow with thermodynamic nonequilibrium between 
the vapor and the liquid. Although research efforts in this aspect is somewhat 
limited because of the difficulty involved, experiments showed that the effects of 
condensation and evaporation cam become of primary importance in determining 
static and dynamic characteristics of saturated vapor bearings and seals. 

The *F-h’ curves for two-phase seal operation obtained by different investi- 
gators has a peculiar feature and is shown schematically in Figure 1.10. The 
positively sloped side of this curve implies a negative film stiffness whereas the 
negatively sloped side means positive film stiffness. For a given dosing force, 
there can be two operating clearances with the smaller one giving rise to un- 





stable operation and the larger one stable. Vaporization also seems to inhibit 
leakage. 

The two models for two-phase seal operation, developed by Beeler and Beatty 
[104,112], work reasonably well at the two extremes - very low leakage rates 
with convection neglected and very high leakage rates with conduction neglected. 
Both models break down as soon as the effect neglected in the respective model 
begins to become important. In actuality, most two-phase seal operations take 
place in the intermediate leakage range when both conduction and convection 
are important. A preliminary model is developed here to bridge the gap between 
the two previous models. This model, known as the ‘Film Coefficient Model’, 
is valid over the entire laminar flow regime unlike the earlier model developed 
by Beeler which only worked at the very low leakage rate end. The new model 
considers both conduction and convection and allows continuous boiling over 
an extended region whereas the earlier model which neglects convection always 
forces a discrete boiling interface and exhibits numerical instability as soon as 
leakage rate starts becoming little higher. With the inclusion of turbulence and 
radial intertia effects, the applicability of the ‘Film Coefficient Model’ can be 
extended to high leakage rate end with the ability to predict choking. Hence this 
model has the potential for describing the seal behavior over the entire range of 
possible leakage rates - low to high. 

Another simplified and semi-analytical model, known as ‘Isothermal Model’, 
has also been developed for low leakage rates. This is based on the model devel- 
oped by Beeler. The assumptions of isothermal condition along the seal interface 
and ideal gas behavior of the vapor permit closed form solutions which may be 
used for preliminary design and analysis. However, to obtain more accurate and 
realistic description, the ‘Film Coefficient Model’ may be used. 

Under certain two-phase operation, seals seem to exhibit self-sustained os- 
cillations even when all the applied conditions remain quite steady. This was 
observed as axisymmetric fluctuations in the film thickness accompanied with pe- 
riodic interface temperature variations. The ‘Film Coefficient’ and the ‘Isother- 
mal’ models are for steady state analysis and none of them can describe this 
anomalous behavior. A preliminary transient heat conduction model has been 
presented here which might help building a comprehensive model to predict the 
oscillatory seal behavior. 

A review of the works done in the single-phase seal area has been presented 
here because some of the ideas and concepts can be used in modeling two-phase 
seal operation. The analysis of two-phase flow in face seals is very complex, and 
a good amount of further research is required to understand a seal behavior. 
The author hopes that the work reported here would give some directions to the 
future efforts. 
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Figure 1.1: Classification Diagram - Industrial Sealing Devices 
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Figure 1.4: Possible Primary-Seal Geometries 
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Figure 1.5: Seal Lubrication Models 
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Figure 1.6: Face Seal with Angular Misalignment and Coning 
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Figure 1.7: Hydrodynamic Patterns on Face Seals 



Figure 1.8: Hydrodynamic Pressure Generation 
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Figure 1.9: Hybrid Gas Seal F-h Curves - Hydrodynamic and Hydrostatic Com- 
ponents 


23 




FORCE 



24 




